First principles calculations of X-ray absorption in an ultrasoft pseudopotentials 
scheme: from a-quartz to high-Tc compounds. 
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We develop a first-principles scheme based on the continued fraction approach and ultrasoft 
pseudopotentials to calculate K-edge X-ray absorption spectra in solids. The method allows for 
calculations of K-edge X-ray absorption spectra in transition metal and rare-earths compounds 
with substantially reduced cutoffs respect to the norm-conserving case. We validate the method by 
calculating Si and O K-edges in a quartz, Cu K-edge in Copper and in La2Cu04. For the case of Si 
and O edges in a quartz and in Copper we obtain a good agreement with experimental data. In the 
Cu K-edge spectra of La2Cu04, a material considered a real challenge for density functional theory 
we attribute all the near-edge and far-edge peaks to single particle excitations. 

PACS numbers: 74.70.Ad, 74.25. Kc, 74.25. Jb, 71. 15. Mb 



I. INTRODUCTION 

With the development of synchrotron radiation 
sources. X-ray absorption spectroscopy (XAS) has be- 
come a very powerful and a widely used technique to in- 
vestigate structural properties and electronic structures 
in condensed matter physics. Since the absorption of X- 
rays at well-suited energies is chemical and orbital selec- 
tive, it is possible to probe electronic excitations and ex- 
plore the local environment around the absorbing atom. 
The use of polarized X-rays allows to separate the con- 
tributions of different atomic orbitals through the study 
of angular dependence of the spectra.^ 

K-edge XAS has been used to study the electronic 
structure of correlated transition metal compoundsf^ii^ii^, 
to probe the local environment around impurities in 
crystals^ and in disordered matter like glass or liquids. 
For instance, XAS plays a crucial role for the understand- 
ing of the microscopic structure of wat er— ili^i^ i ^ ° i ^ ^ . The 
widespread use of XAS as a structural utility and as a 
probe of the electronic structure requires reliable theo- 
retical approaches to interpret the measured spectra. 

Different theoretical methods are available to calcu- 
late XAS. The multiplet approach^^, used to calculate 
pre-edge features for systems with localized final states, 
relies on the solution of a few-sites manybody hamilto- 
nian including several parameters that are suitably cho- 
sen to fit the experimental data. This approach, although 
providing a full manybody solution to the problem, has 
three main shortcomings; (i) it is limited to pre-edge 
structures and (ii) non-local excitations are hardly taken 
into account due to the short-range nature of the consid- 
ered clusters, (iii) in some cases a small variation of the 
hamiltonian parameters leads to substantially different 
spectra. The multiple scattering approachi^illii^ and its 
extension to non muffin-tin potentials have been widely 
used with success, but these methods are not based on 
first principles and require adjustable parameters to in- 
terpret the experimental data. The solution of the Bethe 
Salpeter equation is the first principle method having 



the most satisfactory treatment of many-body effectsi^. 
However being this method extremely time consuming it 
only allows for a description of the pre-edge region of the 
XAS spectra: the near-edge and far-edge regions cannot 
be easily computed. 

Density functional theory (DFT) 

approache o^^'^^i^^'^°i^^'^^ , have been successfully applied 
to K-edges of weakly correlated materials. In a pseu- 
dopotential framework, the use of the PAW^^ method 
allows to reconstruct the all electron wavefunction and 
consequently to obtain XAS intensities unaffected by 
the presence of a pseudopotcntial. Furthermore the 
development of a DFT mcthod^^ using norm-conserving 
pseudopotentials and based on the continued fraction 
approach permits to obtain XAS spectra up to the far 
edge region. Another advantage of DFT approaches is 
that they allow for structural optimization of the local 
environment around the absorbing atom, a key issue in 
the case of impurities or defects. 

In time-independent DFT methods, core-hole effects 
are included in a supercell approach by generating a pseu- 
dopotcntial with a core-hole in the desired atomic core 
level. While this method works very well for weakly cor- 
related system, in the presence of moderate or strong 
correlation it has two main shortcomings. The first is 
the unsatisfactory treatment of electron-electron inter- 
action in the DFT functional. The second is the large 
increase in computational time when dealing with tran- 
sition metals and rare-earths mainly related to the huge 
kinetic energy cutoffs involved and the need to simulate 
large supercells with reduced symmetry .^^ A partial rem- 
edy for the lack of correlation effects is the use of the 
DFT-f U approximation^^. Recently the method of refJi 
was generalized to the DFT-f U approximation-. It was 
shown the DFT-I-U dramatically improves the agreement 
with experimental data in the pre-edge region of corre- 
lated 3d transition metal compounds^. Still the second 
problem holds, namely the huge cutoffs needed to simu- 
late transition metal and rare-earth compounds require 
a substantial computational time. 
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In this work we solve this problem developing a method 
to calculate XAS in an Ultrasoft pseudopotential2i 
scheme and relying on the continued fraction approach. 
US pseudopotentials allow for small cutoffs (20-40 Ry) 
even for transition-metals and rare earth systems, con- 
trary to norm-conserving ones. Thus the use of these 
pseudopotentials reduces the computational cost of the 
supercell calculation by an order of magnitude. The 
drawback is that the continued fraction scheme devel- 
oped in refjii does not apply if at least one ultrasoft 
pseudopotential (not necessary the absorbing atom) is 
present in the calculation. For this reason we reformulate 
completely the continued fraction and the corresponding 
lanczos approach in a way that is suitable for Ultrasoft 
pseudopotentials. We then apply the method to Si and 
O K-edge in a— quartz, to Cu K-edge in copper and in 
La2Cu04 and compare the results with available experi- 
mental data. 

The structure of the paper is the following. In sees. 
ITTl and mil we remind the general expression of the XAS 
cross section within the projector augmented wave for- 
malism. In sec. IIVI and |V] we develop the continued 
fraction approach in the case of ultrasoft pseudopoten- 
tials, and finally in sec. IVII we apply the method to the 
aforementioned systems. 



II. X-RAY ABSORPTION CROSS-SECTION 



1 . 



The XAS cross section isi 



a{uj) = ATT^ahuj2_^\M,^f\^S{Ef 
/ 



E, - hjj) (1) 



where huj is the incident photon energy, a is the fine- 
structure constant, and Mi^f is the transition amplitude 
between the initial state of energy Ei and the final 
state \ipf) of energy Ef. In a single particle approach, the 
many body and jV'/) are replaced by single-particle 
states. Since we consider K and Li edges, l-^i) can be 
either the Is or the 2s atomic core state in the absence 
of a core- hole. The final state \il]f) in the presence of 
a core-hole is obtained in an all-electron first-principles 
calculation. 

In a single particle approach and in the electric 
quadrupole approximation, the transition amplitude is 
given by the matrix element 



with 



l? = e-r-(--(e-r)(k-r) 



(2) 



(3) 



where e and k are the popularization vector and the wave 
vector of the incident beam and r is the electron coordi- 
nate. 



III. X-RAY ABSORPTION CROSS-SECTION IN 
A PAW FORMALISM 

In a first-principles pseudopotential approach, the cal- 
culated wavefunction is namely the pseudowave- 
function of the crystal obtained at the end of the self 
consistent field run. In order to get the all-electron wave- 
functions 1-0/) needed in Eqs. [U [31 all electron recon- 
struction needs to be performed. This is achieved in the 
framework of the PAW method^^ . In this approach the 
all electron wavefunctions are related to the pseudo- 
wavefunctions li/;) though the linear operator T : 



\^) = T\i^) 



(4) 



The T operator is written as a sum of local contribu- 
tions centered around each atomic site R : 



(5) 



R 



The local operators 7r act only within the so-called aug- 
mentation regions centered on atomic sites. Follow- 
ing ref,-i2,, we introduce the all electron (pseudo) partial 
waves I^^R^n) (|0R.n)), and the projector functions (j3r.„| 
that satisfy the conditions'^^ : 



0R,ri(r) 
PR,n 0R', 



i>-R.,n{r) outside Or 



3RR 



(6) 
(7) 



The wavefunctions |(/)r,„) and \4>-R.,n) respectively form a 
basis for valence states, which means that any function 
Xr that vanishes outside is expanded as : 



XR 

Then the operator T is written as 

r = i + ^(l0R,„) 

R,n 



Ixr) 



R,i 



(8) 



(9) 



Substituting eq. [9] in eq. [4] and eq. |4] in eq. [2l leads to 



R,i 



(V'/ PR,n) (^R 



R,i 



V 



(10) 



Since the initial wavefunction ipi is localized on the ab- 
sorbing atom (located at Ro), the terms having R ^ Rq 
can be neglected in eq. [10] to obtain : 



»Ro 



with 



Ro) = X! IP'-o^«) ('^ro,nP|V'i) 



(11) 



(12) 
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Replacing this matrix element in the XAS cross-section 
leads to: 



5{Ef - E, - hw) (13) 



Thus Eq. [T3] express the XAS cross-section in terms 
of single particle states obtained from a pseudopotential 
calculation. Note that in Eq. [12] there is an infinite 
number of projectors. Practically only a few projectors 
are needed to achieve convergence. 



IV. XAS IN AN ULTRASOFT 
PSEUDOPOTENTIAL SCHEME 

In an ultrasoft scheme the norm of the pseudo partial 
waves are different from the norm of the corresponding 
all-electron partial waves. For this reason it is customary 
to define^^ the integrated augmentation charges (ZR^nm 
as: 



iR,r! |</>R,m) " ('/>R,n |0R,m) 



(14) 



The the S operator defined in the ultrasoft scheme^i is 
then: 

5* = 1 + ^ |pRri)'?R,rim(PRm| = 1 + ^ Qr (15) 



R,m,n 



R 



The pseudo hamiltonian H and the pseudo eigenfunc- 
tions 1-0/} satisfy the following equation— : 



H 



ijjf 



EfS 



ijjf 



(16) 



Multiphcation of Eq. [T6]by S^^/^ jj^^ds to 



5-1/2^5-1/2^1/2 



^f)^EfS^'^\i^f) (17) 



The following identity holds (for a proof see app. E]) : 

= lim3[G(x)] (18) 



f 

where x is a real number and 



G(x) - 5-1/2 



-5-1/2 



(19) 



Using Eqlin] and [THl the XAS cross section can finally 
written in a suitable form for a standard Lanczos proce- 
dure : 



cr(Lo) — Anahu! lim S 

7^0 



»Ro 



G{lj + E, 



ORc 



(20) 



where Ei is the energy of the initial state that in a 
pseudopotential scheme is undetermined up to an overall 



constant. In the case of a unit cell having multiple ab- 
sorbing sites which are equivalent under the point group 
symmetry of the crystal, Ei is the same for all the ab- 
sorbing atoms and the choice of Ei simply corresponds to 
a rigid shift of the overall spectrum. On the contrary, in 
the case of nonequivalent absorbing sites in the unit cell, 
the value of Ei depends on the absorbing site due to the 
core-level shift. In this case the choice of Ei is not arbi- 
trary and a careful determination of the core-level shift 
is neededS^. For simplicity in this work we consider only 
examples in which there are only equivalent absorbing 
sites in the unit cell. The determination of the core-level 
shift in the case of multiple nonequivalent absorbing sites 
will be given elsewhere. Thus in what follows we choose 
the energy Ei to be the Fermi level, in the metallic case, 
the highest occupied state, in the insulating case. 



V. LANCZOS PROCEDURE 

Eq. [20] can be calculated using the the Lanczos re- 
cursion method21i2^iSai^. The quantity (0RjG(£'i + 
huj)\(j)-R„) is evaluated using the continued fraction : 



'Ro 



GiE) 



SRc 



»Ro 



»R<: 



ao ^ E ~ i"f 



(21) 



ai —E—if- 



where the real numbers ai and bi are computed recur- 
sively by defining the vectors \ui) such that: 



|mo) 



5-1/2 



"Re 



^Ro 



i'Ro 



S-^'^HS-^'^\u,) = a,\u,) +h,+i\u,+{) +h\u,-i) 
The ai and bi coefficients are defined as: 

ai = (m, 5-i/2#5-i/2 u,^ (22) 
6,, = iu, S-'/^HS-^^^ (23) 



This is essentially a standard lanczos process where the 
initial vector is |uo) and the Hamiltonian H is replaced by 
5-i/2_ff5~i/2. However this is not the more efficient way 
to carry out the lanczos chain since two multiplications 
by 5~i/2 are involved and the S matrix is of the same 
order as the Hamiltonian, namely the dimension is given 
by the number of plane waves in the calculation (the 
kinetic energy cutoff). Thus any application of S~^^^ 
costs as much as the application of H. 

A more efficient way to implement the lanczos process 
is obtained by defining the auxiliary vectors \ti), namely: 



\U) = 51/2 \u,) 



(24) 
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Using this definition, the lanczos process can now directly 
carried out on the jt^) vectors as: 



\to) 



s- 



HS ^\U) ^ a,\U) + b,+i\U+i) +b,\U^i) 

where the new lanczos vectors 1^^} arc no longer orthog- 
onal but {ti\S~^\tj) = Sij. If during the lanczos chain 
the vectors = S^^\ti) are stored then the and bi 
coefficients can be defined as: 



H 



br = (t 



H t 



(25) 
(26) 



Now, each iteration needs only one multiplication by 5'"^ , 
one multiplication by 77, and four lanczos vectors stored 
in memory, namely \ti) and \ti). 

To achieve an efficient implementation of the lanczos 
process, particular care needs to be taken in inverting 
the 5* matrix to calculate 5*"^. Direct inversion of the S 
matrix is unfeasible being the order of the matrix given 
by the number of plane waves. Using the definition of S 
in terms of the Np ultrasoft projectors, the calculation of 
can be performed very efficiently by simple products 
and inversions of matrices of the order of Np x Np , as it 
was demonstrated in refsi^i^. In order to have a com- 
plete description of the method used we recall the main 
passages of the demonstration of ref ^i21 in the appendix 

m 



VI. APPLICATIONS 

The developed method is now applied to silicon and 
oxygen K-edges in a-quartz, to Cu K-edge in copper and 
in La2Cu04. Density functional theory calculations are 
performed using the Quantum-Espresso package^^ and 
the Generalized Gradient Approximation'^^. In the case 
of La2Cu04 we use the Spin Polarized Generalized Gra- 
dient Approximation. The developed continued-fraction 
approach to deal with US pseudopotentials is imple- 
mented in the XSpectra package^^ and distributed with 
the current CVS version of the Quantum-Espresso code. 
Occupied states are eliminated from the spectrum using 
the method of ref,^. The zero of energy is determined 
from the self-consistent calculation on a supercell in the 
presence of a core-hole. In the metallic case we chose the 
Fermi level while in the insulating case the highest oc- 
cupied state. It is important to notice that an insulator 
can become metallic in the supercell calculation due to 
core-hole attraction. In this case the elimination of the 
occupied states is somewhat ill defined, as it is in metallic 
systems. When this occurs, the pre-edge features can be 
incorrect. This is the case in La2Cu04. Further technical 
details of the calculations are given in each subsection. 



A. Si02 (a-quartz) 

Si02 (a-quartz) is a dichroic compound with a hexag- 
onal unit cell and lattice parameters a — 4.9141A and 
c = 5.4060122. The dipolar cross section a has the fol- 
lowing angular dependence^ : 



COS' 



sm' 



(27) 



where e is the polarization vector and 9 is the angle be- 
tween the c axis and e. 

The charge density calculation was performed with a 
2x2x2 supercell containing 72 atoms. Electronic in- 
tegration was performed using only the T point. We 
used a 20 Ry kinetic energy cutoff and a 150 Ry cut- 
off for the charge density, to be compared to the 70 Ry 
kinetic energy cutoff needed in a standard norm conserv- 
ing pseudopotentials calculation^''^ The electronic inte- 
gration in the continued fraction calculation using the 
lanczos method was performed using a centered 3x3x3 
k-points grid of the 72 atoms supercell. Two projectors 
per channel were used in the PAW reconstruction. The 
core-hole width was taken constant and set to 0.8 eV for 
Si K-edge and 1 eV for O K-edge. The continued fraction 
calculation needed around 400 iterations per k-point. We 
performed a calculation with norm-conserving pseudopo- 
tentials in which around 600 iterations were needed for 
comparable accuracy, like in previous worbi. Thus the 
number of iterations needed is smaller, due to the smaller 
cutoff energy. 




-5 5 10 15 20 25 30 35 40 
Energy (eV) 

FIG. 1: Experimental-- and calculated Si K-edge in a-quartz. 
cr|| is the polarization along the c axis, -while ax is the in-plane 
polarization. 

Our results are presented in figures [T] and [H These 
results are in perfect agreement with those presented in 
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FIG. 2: ExperimentaliS and calculated O K-edge in a-quartz. 
(7|| is the polarization along the c axis, while a± is the in-plane 
polarization. 



refiii and obtained with norm-conserving pseudopoten- 
tials. The experimental Si k-edge XAS cross section is 
very well reproduced for both polarizations despite a too 
weak peak around 7 eV. A good agreement between the- 
ory and experiment for O K-edge is obtained. 

The comparison between the ultrasoft pseudopotential 
and the norm conserving pseudopotential calculations on 
a-quartz validates our implementation of XAS using ul- 
trasoft pseudopotentials. 



B. Copper 

Pure copper at room temperature crystallizes in the fee 
structure with lattice parameter 3.6014"^^. The copper 
K-edge XAS cross section was calculated on a converged 
3x3x3 supercell containing 27 atoms. Electronic inte- 
gration was performed over a 10 x 10 x 10 uniform k-point 
grid for both charge density and XAS calculations, with 
a 30 Ry kinetic energy cutoff and a 500 Ry charge-density 
cutoff. 

The calculated XAS cross section is in good agreement 
with experimental data (fig. [3]). The features A, B, 
C, D and E are correctly reproduced, however the peak 
A (shoulder) is shifted of about 1 eV to higher ener- 
gies. This shift is also present in finite difference method 
calculations^. Since copper is metallic, the description 
of the A peak is difficult, because the core-hole attraction 
can drag some states below the Fermi energy. 
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FIG. 3: Calculated Cu K-edge in copper compared to ex- 
perimental data from ref.''-. The Lorentzian 7 broadening 
parameters is variable and runs from 1 to 4 eV. 



C. La2Cu04 

La2Cu04 is the parent compound of high Tc supercon- 
ductors. It is an antiferromagnetic correlated insulator 
considered a challenge for density functional theory. Fur- 
thermore it requires large cutoff energies to be simulated 
with norm-conserving pseudopotentials. Thus it is an 
ideal test for our approach. 

At low temperatures La2Cu04 present an weak or- 
thorhombic distortion of the tetragonal structure. In our 
calculation we neglect the orthorhombic distortion and 
consider the tetragonal structure having a = 5.357 Aand 
c = 13.143 A^. Under this assumption, the Cu K-edge 
XAS dipolar cross section can be described as a linear 
combination of the cross-section having in-plane polar- 
ization {(T±) and of that having polarization along the 
c-axis ((Til), as expressed in eq. [27l 

We treat correlation effects the framework of the the 
spin polarized GGA-fU^i^ approximation, where U is 
the Hubbard parameter on Cu 3d states. The U param- 
eter is 9.6 eV, as calculated from first principles using a 
linear response scheme^i^. We use ultrasoft pseudopo- 
tentials for all atomic species leading to a 30 Ry kinetic 
energy cutoff and a 200 Ry charge density cutoff. The ki- 
netic energy cutoff used with ultrasoft pseudopotentials 
has to be compared with the more than 150 Ry kinetic 
energy cutoff needed in the case of norm conserving pseu- 
dopotentials. We use two PAW projectors per channel 
and non-linear core correction in the Cu pseudopoten- 
tial. By calculating XAS on the antiferromagnetic unit 
cell we have checked that the inclusion of semicore states 
does not affect the result. We used a 1 x 1 x 1 super- 
cell of the antiferromagnetic crystal cell containing 14 
atoms. For the electronic integration we use a uniform 
6x6x6 k-point mesh both for the charge density and for 
the continued fraction calculation. We have verified that 
the result is unaffected by the use of larger supercells. 
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Experimentally, La2Cu04 is an insulator with a gap 
around 2 eV— , and exhibits an antiferromagnetic order— 
with a magnetic momentum on copper atoms around 
0.5^B- Our CGA+U electronic structure calculation 
gives a 0.5 eV charge-transfer gap and a magnetization 
of 0.58^5. 
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FIG. 4: Experimental^ and calculated Cu K-edge in 
La2Cu04. The calculated cross section was obtained with 
a U=9.6 eV parameter on Cu 3d states. The Lorentzian 7 
broadening parameters varies linearly between 1 eV in the 
pre-edge region and 4 eV in the far-edge region. CT|| indicates 
the polarization along the c axis while a± the in-plane one. 
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The results of the CGA-}-U Cu K-edge XAS calcula- 
tions are presented in fig. (3) When the polarization is 
parallel to the Cu02 planes the energy position of the 
different peaks is well reproduced by our calculation. In 
particular B, D,and F are at the correct energy position, 
however the intensity of peak D is underestimated. In 
the pre-edge region the peak A is shifted to higher en- 
ergy. The A peak is not well described in our calculation 
since, due to the underestimation of the electronic gap 
and overestimation of core-hole attraction, the system 
becomes metallic when a core-hole is included in the cal- 
culation. In particular, while the system in the absence 
of a core-hole is insulating, when we add the core-hole 
in the supercell calculation, we obtain a metallic system. 
As a consequence it becomes impossible to distinguish 
between occupied and empty states. The intensity of 
peak A is crucially affected. The system is formally in a 
|3c?^"L) state with no empty d-states and consequently no 
quadrupolar pre-edge. On the contrary it is known that 
a weak quadrupolar pre-edge is present in experiments^. 

The nature of peak B has been widely discussed. It has 
been alternatively assigned to shakedown ls^3(PL — > 
Is^Sd^^L proccsscK^^, and to empty p states of the ab- 
sorbing atom*^-'*'^. Here we unambiguously attribute the 



FIG. 5: Angular dependence of La2Cu04 Cu K-edge XAS 
dipolar cross section compared with experimental datck^. The 
angle between the CUO2 layer and the c-axis is labeled 6, 
so that 9 = 90° corresponds to e along the c-axis. The 
Lorentzian 7 broadening parameters is 0.8 eV. The core-hole 
width has been artificially reduced in the calculation to show 
the presence of different peaks. 



peak B to empty p states and thus this excitation is single 
particle in nature. 

In the spectra having polarization along the c-axis 
the intensities are in better agreement with experimen- 
tal data. However the peak E is substantially shifted to 
higher energies (3 eV). More insight on this issue can be 
obtained by considering the angular dependence of the 
spectra and by comparing it to available experimental 
data^, as shown in FigO 

In the theoretical calculation of Fig. [S] we have sub- 
stantially reduced the peak linewidth to 0.8 eV to em- 
phasize the multi structured form of the different peaks. 
In particular it is seen that the peak at 20 eV (labeled 
peak E in Fig. [5]) is actually composed of two different 
peaks. The high energy one is at the correct energy posi- 



7 



tion while the low energy one is shifted to higher energy 
with respect to the experiment. The incorrect position 
of the lowest energy peak substantially affects the inten- 
sity in the 10-20 eV region for both polarizations and is 
responsible for the disagreement with the experimental 
data. Indeed if this peak was at low energy the intensity 
of peak D in Fig. 2]) would also increase and a better 
agreement with experiment would be obtained for both 
polarizations. 

Kosugi et al4& attributed the decoupling of pre-edge 
and near-edge features present in K-edge XAS on powder 
samples to charge-transfer multi-determinant effects. In 
particular it was shown that the main edge and the B 
peak are decoupled in the isotropic spectra. However the 
pre-edge and edge structures occur at different energies 
when the polarization is in the plane or out of plane. 
We find that the decoupling indeed originates from the 
average over different polarization of the single-particle 
spectra. Thus the doubling of edge and pre-edge peaks 
are not due to multi-determinant charge transfer effects. 

Tolentino et al.~ assigned peaks C and E to 3d}'~'L 
transitions, and peaks D and F to 3cPL transi- 
tions, implying that the XAS spectrum includes multi- 
determinant effects. If this would be the case, being our 
calculation single determinant in nature, a single peak 
should be seen in both directions. This is however not the 
case and our calculation correctly reproduces the main 
experimental features except for an energy shift for the 
peak E. Consequently the near-edge and far-edge struc- 
tures detected in K-edge XAS of La2Cu04 are all single 
particle in origin. Concerning the E peak, it corresponds 
to a single particle excitation that is 3 eV shifted with 
respect to experiments. This shift can be due to an in- 
correct description of the hybridization between Cu 4p 
states of the absorbing atom and La states. La states are 
indeed hard to describe in a single particle approach for 
their intrinsic correlated nature. We believe that, despite 
the 3 eV shift of the E-peak, the single particle origin of 
the C,D,E,F peaks is definitely clarified. To have better 
insight on the subject it would be interesting to study the 
case of Ca2-2;Cu02Cl2^^ since no rare-earths are present 
in the system and consequently the E peak should be at 
the right position. 



VII. CONCLUSION 

A DFT-based continued fraction method using ul- 
trasoft pseudopotentials to calculate the X-ray absorp- 
tion spectra is presented. Our implementation, relying 
on ultrasoft pseudopotentials, is an order of magnitude 
faster than preceding implementations based on norm- 
conserving one^iii. Indeed the bulk of the calculation 
is the determination of the self-consistent charge density 
for a supercell in the presence of a core-hole. Since Ul- 
trasoft pseudopotentials allow for substantially smaller 
cutoffs, the computation cost is strongly reduced. Fur- 
thermore the number of iteration in the continued frac- 



tion is reduced and convergence is then faster than the 
norm-conserving case. 

We validate the method by calculating silicon and oxy- 
gen K-edges of alpha-quartz, Cu k-cdgc in bulk metallic 
copper and Cu K-edge in La2Cu04. In the case of weak 
to intermediate correlation (silicon and oxygen K-edges 
of alpha-quartz and Cu K-edge in metallic copper) we 
obtain a good agreement with experimental data. The 
description of XAS spectra of strongly correlated com- 
pounds as La2Cu04 is typically considered a challenge 
for DFT-based method. Nevertheless we were able to at- 
tribute all the single particle peaks (B, C, D, F in Fig. 
|4]). We then solve the long-standing^i^i^ discussion on 
the attribution of the near-edge and far-edge features in 
La2Cu04. 
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APPENDIX A: PROOF OF EQ. [H 

We proof that for x real number the following holds: 



lim 3 



1 



-5-1/2 



17 



(Al) 



Since 



^Y\i'f)^{Ef -x) ii'f 



7^0 



X — Ef — i'-f 



(A2) 



we have 



1 



/2 



X — Ef — 
1 



X — Ef — i'-f 
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5I/2 



1 



51/2 



y 5-1/2 ^ 



5-2- 



1 



X — S 2 HS 2 — 



^"^E^I^/)(^/ 



= 5-1/2- 



1 



5-i/2i75-i/2 - 



-5-1/2 
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property^: 



where in the last equality we used the following 

— , , . , ^^^^ 
f 

Eq. I All follows from eq. 



APPENDIX B: CALCULATING THE S'^ 
MATRIX 

Following ref.^**, the S matrix can be written as : 



S = ^ + ^qi] \Pi) {Pj 



(Bl) 



The S ^ matrix satisfies the equation SS -'^ = 11: 



SS-^ = (l + ^(7y|p,)(pj|)(Il + ^a,,^|pz)(p„|) 



a + E Ip') ^PjI ('^y + + E ivP^aim) 

(B3) 



Im 



where Pji = {pj\pi)- In matrix form the equation is: 



q + a + qPa ^ (B4) 



where i and j are cumulative indexes for Rm and Rn, 
respectively. We assume that can be written as: 



(B2) 



whose solution is a = —(1 + qP)~^Q. Thus can be 
calculated by inverting matrices of the size of Np x Np, 
where Np is the number of ultrasoft projectors. A similar 
procedure was used in ref.'^^''^^. 
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